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ABSTRACT 

NOAA Active Region (AR) 11429 was the source of twin super-fast Coronal Mass Ejections 
(CMEs). The CMEs took place within a hour from each other, with the onset of the first taking 
place in the beginning of March 7, 2012. This AR fulfills all the requirements for a “super 
active region”; namely, Hale’s law incompatibility and a (5-spot magnetic configuration. One 
of the biggest storms of Solar Cycle 24 to date {Dst = — 143nT) was associated with one of 
these events. Magnetic Flux Ropes (MFRs) are twisted magnetic structures in the corona, best 
seen in ^10 MK hot plasma emission and are often considered the core of erupting structures. 
However, their “dormant” existence in the solar atmosphere (i.e. prior to eruptions), is an open 
question. Aided by multi-wavelength observations (SDO/HMI/AIA and STEREO EUVI B) and 
a Non-Linear Force-Free (NLFFF) model for the coronal magnetic field, our work uncovers two 
separate, weakly-twisted magnetic flux systems which suggest the existence of pre-eruption MFRs 
that eventually became the seeds of the two CMEs. The MFRs could have been formed during 
confined (i.e. not leading to major CMEs) flaring and sub-flaring events which took place the 
day before the two CMEs in the host AR 11429. 

Subject headings: Sun: General — Sun: Surface Magnetism — Sun: Coronal Mass Ejections — Sun: 
Activity — Sun: Corona — Sun: Photosphere 

1. Introduction 

Coronal Mass Ejections (CMEs) are the most 
energetic phenomena of the solar corona. They are 
explosive ejections of coronal plasma that advect 
coronal magnetic flux into the heliosphere, which 
then drives space weather at Earth. We know sev¬ 


eral things about the structure of the advected 
magnetic field. About 40% of CMEs exhibit the 
structure of a Magnetic Flux-Rope (MFR or flux- 
rope thereafter) as determined from coronagraphic 


observations close to the Sun (Vourlidas et al. 


20131. In-situ 1 AU observations suggest that at 


least a third of CMEs maintain an overall MFR 
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morphology as they propagate into interplanetary 
space (e.g. Kilpua et al.||2013 ). Both percentages 
are lower limits. The MFR morphology can be 
easily obscured by projection effects in imaging ob¬ 
servations. In-situ MFR determinations, or Mag¬ 
netic Clouds (MCs) in in-situ parlance, are even 
trickier due to the local nature of in-situ measure¬ 
ment. Several studies argue that the MCs criteria 
are too strict and/or the in-situ detectors inter¬ 
cept CMEs far from their MFR core, say close to 
the CME legs ( Kilpua et al.|2013 Demoulin et al. 
20131. In addition, Riley & Richardson (2013) 


have shown that the propensity of a CME to show 
MFR characteristics decreases with increasing he¬ 
liocentric distance. MFRs are magnetic structures 
comprised of a family of magnetic field lines wrap¬ 
ping around a central field line or axis (hence the 
term flux-“rope”). The fact that MFRs exhibit a 
twisted magnetic field line topology is indicative of 
available free energy and magnetic helicity stored 
within these stressed magnetic structures. 

There is a continuing controversy on the subject 
of CME initiation. The theories suggest that the 
core structure of the CME, the MFR, is formed 
either on-the-ffy (e.g.. Lynch et al.|2008) or exists 


before the eruption (e.g., Kliem & Torok 2006). 


The debate translates into whether the physical 
process behind the eruption is an ideal MHD pro¬ 
cess (pre-existing MFR), or a non-ideal/resistive 
process (on-the-fly MFR formation). Further com¬ 
plications arise from the variety of different for¬ 
mation mechanisms of MFRs. These include flux 
emergence, flux cancellation, magnetic reconnec- 
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2013). Moreover, the pre-eruptive structure may 


be formed at different layers of the solar atmo¬ 
sphere, ranging from the photosphere/lower chro¬ 
mosphere to the corona, or even bodily emerge 
from the convection zone. 

Since the beginning of the SDO era, observa¬ 
tional analyses of pre-eruptive configurations are 
typically case-studies of MFR “candidates” seen 


close to, or at, the solar limb. Cheng et al. (2011) 


have observed indications for the formation of a 
MFR during the impulsive phase of a CME with a 
multi-thermal nature, possibly due to new poloidal 


field added around the MFR via reconnection with 
the ambient field. Zhang et al. (2012) provided 


the first pre-existing MFR observation close to the 
time (a few minutes) of a CME onset. In this ob¬ 
servation a coiled structure, or “hot channel” (de¬ 
tected in 94 A & 131A EUV emission with plasma 
temperatures of « 6.4 & 10 MK) was seen to pro¬ 
gressively unwind into a semi-circular shape and 
then erupt as the core of a CME. This observa¬ 
tion suggested that an ideal MHD instability, like 
the Torus Instability, is responsible for the CME 


eruption. Cheng et al. (2013) determined that the 


“hot channel” was moving faster than the leading 
front of a CME, at least in the beginning of the 
eruption, providing evidence for the role of MERs 


as a driver of the CME. Patsourakos et al. (2013) 


have provided the first direct evidence of a truly 
pre-existing MER. In their observation, the MFR 
was seen after a confined flare, at « 7 hours be¬ 
fore the CME. They proposed a scenario based 
on MFR formation during confined flaring events. 
A recent statistical study using AIA data showed 
that hot MFRs are a common occurrence during 


M and X-class flares, confined or eruptive (Nindos 
et al.|[2M5 ). 


Soft X-ray (SXR) observations can also sup¬ 
ply important information for MFRs within ARs. 


Green et al. (2011) found evidence of flux rope 


formation by means of photospheric cancellation 
in an active region in association with the forma¬ 
tion of sigmoids seen in the Soft X-Rays (SXR) 
and in the EUV light. The sigmoidal structures 
are features which are often seen in SXR prior to 


the eruption (( 

lanfield et al. 

1999 

McKenzie & 

Ganfield|2008 f 

lavcheva & van Ballegooijen|2009 1. 

Rust & Kumar 

(1996) studied a large sample of 


bright sigmoid structures seen in SXR observa¬ 
tions by Yohkoh observatory and they were the 
first to suggest the possibility that SXR sigmoids 
might signify the presence of MFRs. [Green fc| 


Kliem (20091, studying the SXR sigmoid topology 


of an AR, found traces of a bald patch separatrix 
surface which they interpreted as evidence for the 
existence of MFR prior to the onset of a CME. 

The increasing number of reported MER “can¬ 
didates” refer to coronal observations away from 
disk center. The lower background EUV emis¬ 
sion, i.e. there are fewer bright loop-like structures 
along the LOS when MFRs are observed closer to 
the limb, facilitates the detection of the optically- 
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thin MFRs. But to uncover the magnetic struc¬ 
ture, good knowledge of the photospheric mag¬ 
netic field configuration, is necessary. This infor¬ 
mation allows to, at least, pinpoint the “roots” 
of these line-tied 3D MFR structures. In large 
viewing angles, photospheric magnetic field mea¬ 
surements suffer from projection effects (and in¬ 
creased noise) and yield no reliable information 
for ARs very close to the limb. Hence, the past 
studies of the dynamics and formation of MFRs 
were quite restricted because they could not make 
use of coronal magnetic field models, such as the 


Non-Linear Force-Free Field (NLFF, e.g. Canon 


et al. 

2009 Wiegelmann & Sakurai 2012) 

or the 

flux-rope insertion method (e.g. van Ballegooijen 

2004 

Savcheva & van Ballegooijen|2009t Su et al. 

2009 

Savcheva et al.|2012 Bobra et al.|200 

B). For 

instance, in the work of Patsourakos et al. 

20131, 


the event was at the limb and therefore it was 
impossible to support the “confined flare-to-MFR 
creation” formation scenario with magnetic field 
extrapolations. 

In fact, there have been some studies that 
used NLFFF extrapolations to infer the role of 
ideal MHD instabilities of pre-eruptive structures. 


For instance, 

Guo et al. 

( 

2010 

) found indica- 

tions for Kink Instability 

n 

a confined eruption. 

Also, 

Inoue et al. 

(2011 

2012 

20131 analyzed 


the pre-eruptive configuration from NLFFF ex¬ 
trapolations by characterizing the twist of individ¬ 
ual fieldlines numerically, based on the theoretical 


work of Berger & Prior (2006). They found indi¬ 


cations for the Torus Instability as the trigger of 
eruptive events. 

Here, we present unique SDO observations and 
magnetic field analysis that suggest the forma¬ 
tion of MFRs at least 12 hours before they erupt 
as CMEs. In other words, confined flaring/sub¬ 
flaring builds up stressed magnetic fields, which 
eventually give rise to MFRs that erupt. Pat-| 


sourakos et al. (2013) and Tziotziou et al. (2013) 


proposed this confined-flare-to-MFR formation 
scenario solely from observations without support 
from NLFFF extrapolations. 

According to the confined-flare-to-MFR forma¬ 
tion scenario, reconnection takes place first in a 
non-eruptive manner; this changes the topology, 
i.e. arcade to flux rope conversion, and heats the 
plasma within the flux rope structure to flare tem¬ 
peratures. This is why we can observe the hot 


flux ropes during confined flares. Thus the flaring 
serves as a marker of topology change and illu¬ 
minates the resulting magnetic structures. Flux 
cancellation as described in [van Ballegooijen &:| 


Martens 

( 

1989 

) or 

Green et al. 

( 


also 


volves magnetic reconnection and arcade-to-flux 
rope conversion, but requires the disappearance of 
magnetic elements in the photosphere; a process 
which is not always observed. 

The source region was d-spot NOAA AR11429 
(for a full-disk AIA context image see Figure]^. 
AR11429 was a “super-active” AR, which ejected 
several CMEs and a multitude of flares during its 
transit over the solar disk. We focus on two fast 
CMEs (> 2 000km/s) which occurred within one 
hour from each other on March 7, 2012. These 
CMEs, which were initiated from different parts of 
AR11429 and followed different propagation paths 
in the low corona and beyond, are of general inter¬ 
est as they were associated with one of the most 
intense geomagnetic storms of cycle 24 to date 
{Dst = —143nT; Richards^|2013 ). 

The paper is structured as follows; in section 2 
we provide an overview of the observations. In 
section 3 we discuss the coronal activities from 
the time of the confined flaring to the time of the 
twin CME eruptions. The NLFFF extrapolation 
specifics and the associated analysis and findings 
from magnetic data products are discussed in sec¬ 
tion 4. In section 5 we discuss the evidence sup¬ 
porting the “dormant” MFR scenario. 

2. Overview of Observations 


We use high temporal cadence and spatial res¬ 
olution observations taken in different spectral 
bands by the Atmospheric Imaging Assembly 
(AIA; [Lemen et al. 2011) and the Helioseismic 
and Magnetic Imager (HMI; Scherrer et al.|[2M2 ) 
aboard the Solar Dynamics Observatory {SDO', 
SDO orbits Earth at an 
For three weeks 


Pesnell et al. 2012). 


inclined geosynchronous orbit, 
around the Vernal/Autumnal equinox, the so¬ 
lar disk is briefly obscured by the Earth on a 
daily basis, in what is known as the SDO eclipse 
season. Since the events under study were in 
early March this introduced daily data-gaps in 
our observations. For a different vantage point, 
we used images from the Extreme Ultraviolet Im¬ 
ager (EUVI B for short; Wuelser et al.||2004 | at A 
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195 A part of the Sun-Earth Connection Coronal 
and Heliospheric Investigation instrument (SEC- 
CHI; [Howard et al.|[2008| aboard the Solar Ter¬ 
restrial Relations Observatory (STEREO; Kaiser 
et al. 20081. The STEREO mission comprises 


two satellites, one moving ahead (STEREO “A”) 
and one behind (STEREO “B”) the Earth. Only 
STEREO-B was appropriately positioned at that 
time. At the time of the observations STEREO 
B was at « 117° east of the Sun-Earth line. Both 
AIA and EUVI-B as well as HMI are full-disk im¬ 
agers of the coronal plasma and the photospheric 
magnetic field respectively. 

2.1. Magnetic Field Observations in the 
Photosphere 

Eor the vector magnetic field observations we 
used preprocessed cutouts (containing the AR) at 
full-cadence (720 sec). The preprocessing consists 
of the resolution of the 180°-ambiguity for the az¬ 
imuthal component of the magnetic field, and the 
transformation of the helioprojective images to a 
Cylindrical Equal Area (CEA) projection, which 
preserves the pixel size (Hoeksema et al. 20141. 


The CEA projection is essentially a transforma¬ 
tion of the observational data into a Cartesian 
geometry, as if we were observing the AR “from 
above” while commoving with its guiding center. 
The resulting CEA dataset comprises three sets 
of 2-D Cartesian images, that is, one image for 
each component of the photospheric vector field, 
Bphot = (BP^°*, Bf^°*). The physical pixel 

size in the CEA dataset is 0.36Mmpixel“^ and 
the size of the original cutouts was reduced to 
205Mmxl45Mm (or 564x391 pixels) which en¬ 
veloped the AR. Due to the eclipse season, the 
total number of vector cutouts is reduced to 228 
frames (instead of 240). A sample vector image 
from this time-series (showing all three compo¬ 
nents of Bphot) is presented in Figure]^ 

Figure is illustrating the complexity of 
AR11429. It is an anti-Hale AR, which means 
that the leading polarity is opposite to that of 
the hemispheric trend. In fact AR11429 is the 
strongest (in terms of total unsigned magnetic 
flux; $ = 5.9 X 10^^ Mx) anti-Hale AR of Cycle 
24 to date. It exhibits a complex, d-spot con¬ 
figuration (specifically, it is /J/y/d-configuration). 
<5-spots are known for their extreme activity (e.g. 
jZirin fc Tanaka|[l973[ |Tanaka||l991| |Fan||2009| and 


references therein). 

Another interesting feature of AR 11429 is the 
existence of a complex, sharply-defined and evolv¬ 
ing PIL as seen in the Bz component of the pho¬ 
tospheric magnetic field. In particular, the west¬ 
ern part of the AR exhibits multiple PILs (Fig[^ 
SW Box), in contrast to the eastern part of the 
AR which features a simpler and longer PIL. So 
we have the rare occurrence of PILs of different 
spatial scales coexisting within the same AR, and 
all of them exhibiting strong spatial gradients (i.e. 
very dense PILs). 

A useful metric for the strength of the spatial 
gradient associated with PILs is the R-parameter 
(Schrijver 20071. It parameterizes the unsigned 


magnetic flux near high gradient, strong field PILs 
from LOS magnetograms. In our case, we get a 
logioR ~ 5.2 which lies at the high-R tail of Schri- 
jver’s statistical study. However, the exact value of 
R might differ depending on the instrument (MDI 
vs HMI) and their properties (filling factor, resolv¬ 
ing power and point spread function). Neverthe¬ 
less, AR 11429 has a strong gradient PIL which is 
quite often seen in d-spot ARs. Also, in Figure 
one can see that the horizontal field vectors (shown 
with blue color) close to the PIL are oriented along 
the PIL. This is evidence of strong non-potential 
shear and we discuss this in more detail in Section 

imi 

In Figurej^we present the time evolution of the 
photospheric vector field during March 6, 2012. 
It is quite obvious that significant changes occur 
along the PIL by means of shearing motions of in¬ 
dividual flux elements (green circles; most notably 
in the NE side of the PIL - white “blobs”). These 
“blobs” originate from a small emergence event 
on March 5, 2012, 10:30 UT close to the Northern 
negative polarity of the AR. In turn, the shape of 
the PIL (shown in red; found for every frame by 
means of an image gradient operation as injZhan^ 


d(p 


et al. 2010) is changing as this can easily be seen in 


Fig[^ Shearing also occurs at the location of the 
SW negative sunspot with the elongated positive 
polarity on top if it. 

Another interesting trend of the photospheric 
evolution is the convergence of the southmost neg¬ 
ative sunspot towards the PIL. In addition to all 
these motions, AR 11429 displays sunspot rota¬ 
tion, which cannot be easily appreciated in Fig|^ 
To better appreciate the reported PIL dynamics. 
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a movie of CEA vector frames spanning a period 
of 5 days is available in the online version of the 
paper (filename “moviel.mpg”). To summarize, 
these observations suggest that the shearing mo¬ 
tions and sunspot rotation contribute to the build¬ 
up of the non-potential shear angle in the photo¬ 
sphere, and of stressed magnetic held conhgura- 
tions in the corona. The importance of these ob¬ 
servations will become apparent in Section 4. 


2.2. Magnetic Field Extrapolation 


Due to the on-disk location of AR 11429, it is 
possible to use photospheric magnetic held extrap¬ 
olations to assist our investigation on the structure 
of the coronal magnetic held. A reasonable ap¬ 
proximation for the magnetic helds in the corona 
is the Force-Free assumption, i.e. the absence of 
Lorentz Forces by neglecting the cross-held electric 


currents (for a discussion see review by Wiegel- 
mann fc Sakurai|2012 ). In other words, the model 
requires that the electrical currents how strictly 
along the magnetic held lines. Another require¬ 
ment is the absence of magnetic monopoles in the 
bounded volume, V. Using vector notation, the 
above can be expressed as 


V X B = a{x,y)B ; V • B = 0 (1) 


where a{x,y) is the force-free or torsion parame¬ 
ter, which in general is a function of position but 
is conserved along each held line. The above set of 
equations is describing the general case, i.e. that 
of Non-Linear Force-Free Fields (NLFFF), which 
is considered quite successful in reproducing the 
overall morphology of ARs and current channels 
along PILs ( Wiegelmann fc Sakurai|[2012 and ref¬ 
erences therein). A special case is the linear force- 
free case, when a(x, y) = const globally, which 
is known to be inaccurate when the extrapolated 
magnetic held lines are compared to structures 
from EUV images ( Wiegelmann fc Sakurai|[2012 1. 

For our investigations we used an optimization 
technique (see Wiegelmann 2004) for computing 
the NLFFF in the corona, based on the method 
introduced by Wheatland et al. (20001. Using ap¬ 
propriate boundary conditions (see following sub- 
and Appendix A), this numerical 


2 . 2.1 


section 

method yields a NLFFF solution by minimizing 
a penalty function, L, in the computational vol¬ 


ume, V, as 


L = j w(x,i/,z)[R-2|(VxB)xB|2-p|V-Bp]dU 
y 

( 2 ) 

where w{x^ y, z) is a scalar function with a value of 
1 in the physical domain of the volume and drops 
smoothly to zero when approaching the top and 
lateral boundaries. Obviously, when the penalty 
function L = 0, both the Lorentz force is zero and 
the solenoidal condition is satished in the entire 
computational volume, V, which now contains the 
NLFFF. 

2.2.1. Input Data to the NLFFF Code 

The input data for the NLFFF calculation were 
photospheric vector magnetogram cutouts from 
the HMI (discussed in section . We performed 
NLFFF extrapolations at an hourly time-step, 
covering the period from March 6, 2012 00:00 UT 
to March 7, 2012 23:00 UT. Due to the eclipse 
season, extrapolations at 08:00 UT of both March 
6, 2012 and March 7, 2012 were not calculated. 
Moreover, the boundary of March 7, 2012 at 07:00 
UT was determined to be problematic (out-of- 
focus) and it was discarded. We ended up with 45 
NLFFF extrapolations (3-D vector datacubes) per 
time-step. The size of the NLFFF extrapolation 
datacubes used in our analysis was 250x163x106 
pixels (or 180xll7x76Mm^) for each component 
of the NLFF magnetic field (see appendix for more 
information). 

2.3. Coronal Observations 

We analyze the AIA coronal observations from 
March 5, 2012, 23:59:58 UT, when the AR was 
at E43°N18°, to March 6, 2012, 23:59:58 UT 
E17°N18°. The spatial resolution of the AIA im¬ 
ages is 076 pixel”The passbands used (and 
their nominal formation temperature) are: 131A 
(0.4MK & lOMK), 94A (6MK), 211A (2.0 MK), 
I 93 A (1.6MK & 20MK). The time-cadence is 12 
sec, resulting in 14400 images per passband. Two 
data gaps exist (one per day between 06:50 and 
07:40 UT) due to the SDO eclipse season. From 
the full-disk images we extracted a sub-field cen¬ 
tered at the AR at a size of 600" x 500". The Field 
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of View (FOV) size was chosen to fully contain 
the AR, while it rotates at the Carrington rota¬ 
tion rate. 

From the four available passbands of EUVI-B, 
we selected the 195 A passband because it has the 
highest cadence (5 min) and probes the hottest 
plasma. The 195A channel contains a weak Fe 
XXIII line which contributes only during flares 
and is of interest to us since we are searching for 
a hot MFR. The pixel size is 1" 59. There were 
576 images in total and the AR was located close 
to the EUVI-B western limb. The combination of 
SDO/AIA and EUVI-B observations provided a 
comprehensive view of AR 11429 from above and 
from the side, respectively. 

3. Coronal Activity Prior to the CMEs 

As we mentioned before, the transit of AR 
11429 is marked by intense coronal activity with 
many eruptive and confined flare events. During 
the period of interest here, covering the entire day 
before the onset of the CMEs of March 7, 2012, the 
GOES satellite records several SXR flares above C 
level in the 1-8 A X-ray flux (Figure]^. We use 
the 131A images, covering AR 11429, to locate the 
source of this activity. It manifests itself as tran¬ 
sient brightenings, either strong enough to be asso¬ 
ciated with a GOES flare or weaker in the sub-flare 
domain, of hot (^10 MK) coronal structures. The 
close association between the 13lA light-curve of 
AR 11429 and the GOES light-curve suggest that 
the SXR emissions originate from this AR. 

The brightenings occur predominately in two 
locations, in close proximity to the PILs and have 
a very elongated character (Figure]^. To examine 
their relationship to the underlying complex mag¬ 
netic field structure, we formed an AIA/HMI com¬ 
posite movie (filename “movie2.mpg”) using the 
131A and the Line of Sight Magnetogram (BLOS) 
images (blue negative, orange positive). We show 
a set of representative snapshots in Figure and 
the movie is available in the online version of the 
paper. The BLOS background images are satu¬ 
rated at ±30 G in order to serve as a PIL-tracer, 
which is found at the interface between the two 
color-coded polarities; the foreground 131A image 
shows that the transient brightenings reside in the 
intersection of the blue and orange saturated areas 
of BLOS, i.e. in the vicinity of the PIL. The movie 


spans five days and a large range of heliographic 
longitudes. The close proximity of the brighten¬ 
ings with the PIL, over such a long time, is indica¬ 
tive of low-lying active structures. Thus, the 131A 
evolution indicates two main kernels of activity 
within the AR. One kernel is associated with the 
NE PIL and the other with the SW PIL. In fact, 
the SW kernel is active all the time -i.e. bright 
in the 131A passband-, while the NE is active 
but in a more intermittent fashion. In the begin¬ 
ning of March 6, 2012 the NE brightenings seem to 
originate from a fragmented structure, apparently 
small hot-loop-channels tilted with respect to the 
NE PIL. The structures traced by the transient 
brightenings become more monolithic—they tran¬ 
sition into an elongated, possibly low-lying struc¬ 
ture, almost perfectly aligned with the NE PIL, 
towards the end of March 6, 2012. Another impor¬ 
tant characteristic of the hot-loop-channels seen 
in Figures and 1^ is that the highlighted struc¬ 
tures (either pointed with an arrow or dotted lines) 
have a weak sense of an S-shape and hence maybe 
weakly twisted. 

3.1. Confined Flaring Activity during 
March 6, 2012 

12 hours prior to the GMEs of March 7, 2012, 
an M2.1 flare occurred around 13 UT on March 
6 , 2012 (labeled “confined” in Figure]^. As dis¬ 
cussed earlier, the transient activity, in the form of 
brightenings in 13lA, is continuous and localized 
at the NE and SW PIL (Figures [^ . The M2.1 
flare brightens a large part of the NE PIL and the 
corresponding structure seen in 131A (e.g. top 
right panels of Fig dotted line; Fig red ar¬ 
row) gives the impression that is being lifted up. 
It becomes fainter as it (presumably) cools. The 
same reappears shortly before the onset of the first 
GME of March 7, 2012 and eventually erupts with 
this GME. The location of M2.1 flare lies in the 
middle of the NE PIL, as determined from the 
131A image on March 6, 2012, 12:35 UT and is 
the same location for the eruptive X5.4 flare at 
00:02 UT of March 7, 2012 associated with our 
first GME. 

There were no large-scale EUV dimmings asso¬ 
ciated with the M2.1 flare as can be seen in 211A 
base difference images in Figurej^ A similar result 
was reached with 193A base difference images. No 
GME was recorded in the available coronagraph 
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data (SECCHI, LASCO). The lack of large-scale 
EUV dimmings and coronagraph CME signatures 
lead us to conclude that the M2.1 flare was a con¬ 
fined flare. 


One hour after the peak of the confined M2.1 
flare, a Differential Emission Measure (DEM; fol¬ 
lowing [Plowma^^eF^ 20131 analysis verifies the 
existence of a hot (~10MK) structure in the NE 
PIL (bottom 6 panels of Figure [^. The struc¬ 
ture is also seen in the 131A images taken af¬ 
ter 12:23 UT (Figure [^. We interpret this as 
the formation and filling with hot plasma of an 
MFR structure during the M2.1 flare. The overall 
shape and morphology of the NE structure associ¬ 
ated with CMEl is very similar to the one seen at 
12:23 UT. This is demonstrated in the top right 
and bottom right panels Figure (dotted yellow 
lines are added to envelope the bright structure). 

Detailed spectroscopic diagnostics of AR11429 
by the EIS spectrometer on Hinode, around 12:00 
UT and 21:00 UT on March 6, 2012 also show 
evidence for hot MFR structures (Syntelis, Gon- 
tikakis, Patsourakos, Tsinganos 2015; in prepara¬ 
tion). 

A second M-class flare (Ml.3) occurred around 
21 UT on March 6 (Figure]^ at the same location 
as the first M-class flare. For the sake of brevity, 
we do not describe the analysis details here but 
the coronal response was very similar to the first 
M-class flare - i.e, no CME, a confined flare and 
“illumination” of MFR-like structures. We thus 
have two confined M-class flares in the midst of a 
plethora of transient brightenings resulting in the 
creation or illumination of a hot MFR structure 
at the location of a fast CME 12 hours later. 


Regarding the SW PIL, there is persistent sub- 
flaring activity with small transient brightenings 
which are more spatially constrained than those 
in the NE PIL (due to its smaller lengthscales). 
There were no CMEs originating from that loca¬ 
tion during March 6, suggesting that the observed 
sub-flarings were confined. The activity may indi¬ 
cate the sites of stressed magnetic structures that 
could have contributed in the formation of the 
MFR of CME2 (see Fig|^. This PIL is the lo¬ 
cation of the second X flare on March 7, ~ 50 
mins after the X flare at the NE PIL. 

Observations from the Earth (SDO/AIA) and 
from a different viewpoint (from the 195A chan¬ 


nel of EUVIB) suggest the slow rise of the sys¬ 
tem during the M-class flare of March 6, 2012. To 
measure the speed of the rise we created slit-time 
plots seen in Fig[^ The position and orientation 
of the slit was selected to intercept the images 
along the direction of fastest expanding motions 
(in the online version of the paper a multipanel 
movie is available, “movie3.mpg”). The measured 
speeds of the bright front seen in the right panels of 
Fig|^ are ~ 10-20km/sec (derived from the slope 
of the “step”-like feature in the slit-stack maps 
covering 11:00 - 15:00 UT of March 6, 2012; left 
panels). These speeds represent a small fraction of 
the Alfven speed in the corona, and signify a slow 
quasi-static rise of a magnetic structure possibly 
caused by slow photospheric motions (shearing or 
twist). We essentially have a two-phase kinematic 
behavior around the flare event. First, a quasi¬ 
static slow rise (the “step”-like feature in the slit- 
stack maps) at 10-20 km/s from ~ 12:23 UT dur¬ 
ing the flare. The rise phase ends at 12:40 UT. 
Thereafter we observe much shallower slopes sig¬ 
nifying that the magnetic structure has reached a 
stable position. This is additional evidence for a 
confined flare. 


In the cooling phase of the M2.1 flare (^13:46), 
cool filamentary plasma appears along the pre¬ 
sumed MFR axis. The red arrows in Figure [T0| 
show a dark absorption feature and hence cooler 
than the peak temperature of 195 A, 1.6MK. This 
suggests that a coherent structure exists along the 
PIL to support cool plasma. This observation is in 


agreement with MFR observations of Patsourakos 


et al. (20131 that show cooling at the MFR center. 


Note that the same absorption can be seen from 
the AIA vantage point (almost seen from above, 
top right panel in Figure]^. 


4. MFR Tracers from the NLFFF Extrap¬ 
olations 


To further aid in our search for pre-existing 
MFRs, we turn to the NLFFF extrapolations. We 
use them as a “laboratory” to generate several 
diagnostics to explore the existence of MFRs in 
the pre-eruptive magnetic configuration. In the 
following paragraphs we discuss these diagnostics 
separately. 
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4.1. Magnetic Stream-line and Electric 
Current Visualization 

In Fig we show the extrapolated NLFFF by 
means of magnetic stream lines on March 6, 2012, 
23:48 UT as seen from the top of the domain. The 
time corresponds to shortly before the onset of the 
first CME. For the purpose of clarity and in order 
to give a sense of 3-D depth, we color the magnetic 
stream lines with green color for the overlying field 
and with teal color for the lower lying magnetic 
structures. In the online version of the paper a 
movie showing a fly-by around that extrapolation 
is available (“movie4.avi”). The low-lying mag¬ 
netic structures exhibit MFR morphology when 
observed very close to the PIT and they gradu¬ 
ally become sheared for lines that are tied further 
about the PIT. The overarching, higher lying mag¬ 
netic structure seems to describe well the EUV im¬ 
ages and seems closer to a potential state. As sug¬ 
gested from the 13lA brightenings, the majority 
of the transient activity originates from low-lying 
magnetic structures along the PIL. Strictly speak¬ 
ing, there are two main kernels of activity, one in 
the NE PIL and another in the SW PIL (as seen 
in the previous figures). 

By applying Ampere’s law as V x Bnlfff = 
/roJ numerically, we get the spatial distribution of 
electrical currents J in the NLFFF. Their mag¬ 
nitude, |J|, represents a proxy for ohmic heat¬ 
ing. In other words, regions with strong |J| would 
be more heated and thus exhibit higher tempera¬ 
tures. A more complete approach should incorpo¬ 
rate a full energy equation, which includes source 
(e.g.. Ohmic heating) and cooling (e.g. radiation 
and conduction) terms to work out the detailed 
plasma emission in any AIA channel. However, a 
careful inspection of the 131A image time series 
with images containing the areas of strong electri¬ 
cal currents in the NLFFF datacubes shows that 
the 10 MK plasma emission is described well by 
the location of the stronger NLFFF-inferred elec¬ 
trical currents. The strength of these currents is 
characteristic of the compactness of the PIL. 

We created |J| time series which contain the 
LOS integral of |J| (middle row panel of Fig 121. 


Also we took vertical planar cuts in the middle 
of the NLFFF boxes that show |J| in color con¬ 
tours with the local magnetic field projected onto 
the cut plane (bottom row panel of Fig 12 the 


leftmost “lobe” corresponds to the NE MFR and 
the rightmost to the SW MFR). The oval shape 
of the B-projection suggests the existence of MFR 
structures detached from the photosphere. The 
cores of the structures are low-lying (« 3Mm). 
The length of the two prevalent current channels 
is 80 Mm for the NE MFR structure and 40 Mm for 
the SW. The radius of the MFRs is about 10 Mm 
(NE) and 6Mm (SW) (deduced from where |J| is 
20% of the peak current). The mean total mag¬ 
netic field strength measured on the cut plane of 
the NE MFR (before the eruption) is ~ 650 G and 
for the SW MFR, ~ 500 G (see Table 0. 

There is also a striking simil arit y between the 
131A and |J| time series of Fig[^ At the begin¬ 
ning of our time series (early March 6, 2012) the 
transient brightenings along the NE PIL delineate 
a fragmented active structure. As the time goes by 
and towards late March 6, the NE PIL is under¬ 
going transient brightenings, but now as a more 
monolithic structure. This is evidence for highly- 
sheared (almost parallel to the PIL) structures 
that seem to comprise a very elongated structure 
along the PIL. This could be a result of a continu¬ 
ous cancellation process which creates a very elon- 


gated structure with a strong axial field ( 

van Balle- 

gooijen & Martens||1989 

Green et al.||2011 

). This 


fragmented-to-monolithic current channel forma¬ 
tion is reflected in the NLFFF |J| datacubes. 


Persistent transient brightenings also occur in 
the shorter, SW PIL. There is a very good cor¬ 
respondence of the location of NLFF volumetric 
currents with the location of the brightenings in 
131 A. However the picture is more complicated 
when inspecting the vertical cuts; there seems to 
be a low lying persistent structure, or “island”, 
between the structures associated with the NE 
and SW MFR. This structure exhibits an opposite 
sense of rotation, or chirality, for the magnetic field 
lines. More details can be found in Section lT3l 
The overall picture we get from the location 
of the two strong current channels along the PIL 
and the helical morphology of the magnetic field 
around these current channels is suggestive of the 
existence of two MFR candidates. 













4.2. Twist Number and Mean Twist Evo¬ 
lution 


An alternative way to detect twisted magnetic 
field structures in NLFFF extrapolations was in- 



T„, for magnetic field lines with non-zero helicity, 
as 

(3) 

where, a, is the force-free parameter, and, dl, is 
the infinitesimal length of the field line. For force- 
free helds, a is constant along a field line and the 
twist is then given by 

Tn = ^aL (4) 


where now, L, is the total length along a field line. 
We derive a as V x Bjp/ioi- Thus, we can obtain 
the twist number for each field line since we know 
its length a from NLFFF extrapolation. 


The method by Inoue et al. (2011) is capable 
of producing 2D maps of the twist number, r„, by 
mapping Tn at the footpoints of each magnetic line 
by simply multiplying a-maps and “length-maps”. 
However, caution must be taken since a may not 
be the same at both footpoints. This can be either 
due to instrumental noise, deviations from force- 
free condition or noise due to numerical deriva¬ 
tives associated with the calculation of a. In our 
implementation, we simply take the average value 
of the a between the negative and positive polar¬ 
ities as d = and then we smooth the map 


with a 3x3 gaussian kernel (following Inoue et al. 


20111. For the line integration we use a second or¬ 


der Runge-Kutta integrator with a 1/16-pixel in¬ 
tegration “time-step”. The line integration starts 
from pixels below a threshold < — 500G 

(i.e. from negative footpoints). When the line in¬ 
tegration reaches the positive polarity at the pho¬ 
tosphere, the integration stops and the value of to¬ 
tal length of the field line is stored on a 2-D map 
on both positive and negative footpoints. We call 
this the “length-map”. After the entire negative 
polarity “seedpoints” have been used for the inte¬ 
gration, the map looks smooth in the negative po¬ 


larity but the positive polarity seems “sparse” (i.e. 
contains holes, in other words, points of avoidance 
for the final locations of the integration). This sig- 
nihes that the NLFFF method does not perfectly 
satisfy the solenoidal condition at the boundary, 
thus we resorted to a work-around described in 
the appendix B. 


A major concern using T„ maps for detect¬ 
ing MFR structures is the possibility of confusing 
Sheared Magentic Arcades (SMAs) with MFRs. 
So we need to define T„ thresholds separating 
MFRs and SMAs. For this purpose, we used an 
analytical expression (Aschwanden 2004) to de¬ 
rive various SMA fields and found that they could 
reach a maximum |T„| of 0.5. So we considered as 
an MFR field line any line with |r„| > 0.5. Our 
choice of the absolute twist threshold between a 
SMA and an MFR seems reasonable and conser¬ 
vative. The inspection of the maps ( Fig[^ 
reveals that: (1) we obtain a few locations of 
relatively high twist around and along the PIL 
in accordance with the imaging observations in 
the corona, and (2) \Tn\max is not extremely big 
(not signihcantly exceeding one full turn). These 
findings are indicative of weakly twisted low-lying 
MFRs, in agreement with the EUV observations 
in the hot 131A channel. 


In Figj^we show the temporal evolution of the 
spatial averages of the force-free parameter, twist 
number, and the fraction of twisted flux above se¬ 
lected twist thresholds over total flux, dehned as 


~TBz 


( 5 ) 


The surface integral is taken over the negative po¬ 
larities. 

The analysis is showing an overall increasing 
trend for the average a, T„ and F(r) over the en¬ 
tire considered time period. For the average force- 
free a there seems to be a decreasing trend before 
the onset of the CME events (vertical line at 00:00 
UT) and a rapid increase after the events. The av¬ 
erage T„ shows a steady value of about a quarter 
turn with negative chirality which increases slowly 
after the CMEs. The fractional twisted flux F(r) 
for \Tn\ >0.5 (half a turn) shows qualitatively the 
same behavior to the average twist. However the 
fraction of twisted flux above half a turn is 10% 
and steadily increases after the CMEs. The inter- 
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pretation for this is due the continuous growth of 
the AR in terms of magnetic flux (in the form of 
emergence of small magnetic features throughout 
the time-range of the observations). Also,in ad¬ 
dition to the continuous emergence, the emerged 
flux is (or becomes) more stressed so that to be 
reflected in magnetic held metrics associated with 
its stress (e.g., twist number and alpha parame¬ 
ter). 

Finally, in each of the top two panels of Fig[T^ 
we show the total number of footpoints in the neg¬ 
ative seed- “patch” (which were initially selected as 
< —500 G). There is an increase of the num¬ 
ber of points (from 8700 to 9500 points) after the 
CMEs which implies that there was an increase in 
the magnetic flux at the photospheric boundary 
which supports our interpretations of the previous 
paragraph. 


4.3. Photospheric and Coronal Chirality 
Patterns 

Since the structures undergoin g tr ansient bright- 
enings in 131 A (Figures and 12) are strongly 
emitting and stand out well about the background, 
it is possible to identify the helicity sign, or the 
chirality of the structures and compare it with the 
NLFF results. The chirality can be obtained, as 
a 2-D map, from the sign of the force-free param¬ 
eter, a. For this task, we convert the chirality 
maps from the native heliographic (CEA) coordi¬ 
nates of the NLFFF into helioprojective. Now, we 
can identify where the “active” (in 131A) PIL- 
aligned structures are rooted and thus estimate 
their chirality. Since these are transient brighten- 
ings, their intermittent nature can be summarized 
best in a single map by calculating the standard 
deviation of the 131A time-series for each pixel. 

In Figure [15] we present composites of a-maps 
and 131A standard deviation maps accumulating 
activity for two hours before the confined M-class 
flare and the CME eruptions. Again, the maps 
show two clear clusters of high-temperature coro¬ 
nal activity: one in the NE PIL and the other in 
the SW PIT. The dominant chirality of the NE 
PIL “structure” or cluster in the 131A is negative 
(overlaid on red a-patches) but for the SW PIL the 
picture is complicated. In fact, there is a patch or 
“island” of positive chirality (blue a-patch) in be¬ 
tween the negative chirality of the NE PIL and the 
negative chirality at the SW PIL (see also vertical 


cuts at bottom panels of Fig 12 the positive chi¬ 
rality structure corresponds to the middle current 
structure). 

In addition. Figure shows that this is the 
location of continuous activity of small low-lying 
structures. This holds true throughout our observ¬ 
ing period. Transient brightenings indicate contin¬ 
uous, small-scale reconnection events (no big flares 
during the considered intervals) localized within 
the two distinct regions—one in the NE and the 
other in the SW. Thus, the chirality patterns sup¬ 
port the existence of 2 MFRs in the location where 
the two CMEs occur. 


4-. 3.1. Shear Angle Evolution 

Another way to characterize the deviations 
from potentiality is by means of the shear an¬ 
gle of the photospheric magnetic field on the PIL 


( Hagyard et al.|[l984 Ambastha et al.|[l993 ). The 
shear angle is defined as the angle between the 
observed photospheric field and the potential field 
at the photospheric layer. 


9 = cos 


^o6s * ^ 


pot 


Boh.JIB 


pot I 


( 6 ) 


phot 


The larger the shear angle, the further from 
the potential state the magnetic field is at the 
PIL. The main process responsible for gradual 
shear build up is the slow photospheric motion of 
magnetic elements, which continuously drag (and 
shear) magnetic field lines. However, there have 
been reports for rapid, permanent changes of the 
shear during flares. Their interpretation is unclear 
as both a rise and a drop of the shear angle have 


been observed after the flare (Wang et al. 1994 


Petrie 20121. The horizontal field in the photo¬ 


sphere could become more sheared following the 
formation of a vertical current sheet when a MFR 
moves upward (e.g. Forbes & Priest 19951. In 


our case, we calculated the shear angle 9 using the 
eq (|^ by sampling the shear angle along the PIL. 

The PIL is found for every frame by means of 
an image gradient operation (as in Zhang et al. 


20101 and the end product is a binary map. The 


resulting PIL maps are grown to a thickness of 4 
pixels. A histogram of the shear angle distribution 
versus time is shown in Fig Overall we get 
high shear angles and that is very important for 
the evolution of the AR, since shear can generate 
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twisted field lines. 


In the beginning of March 6, the mean shear 
angle is around 68° and begins to slowly increase 
for about 6 hours before the CME events, reaching 
a maximum of 74° at the moment of the eruption 
(left panel of Fig 16). This number differs from 
the maximum number reported by Petrie| (2012) 
(~ 90°) as they focused on a small part of the NE 
PIL. In the right panel, we show the total number 
of pixels above a threshold of shear angles. It is 
interesting that the number of pixels in the PIL 
with shear angles above 75° is increasing steadily 
until the time of the CMEs; this is understood as 
there are more magnetic vectors in the PIL ori¬ 
enting along the PIL since if the structure was po¬ 
tential, the magnetic vectors would be vertically 
oriented to the PIL. These measurements indicate 
the presence of magnetically stressed structures in 
the locations of the two CMEs, which further sug¬ 
gests the existence of 2 parent MFRs. 


4-3.2. Critical Decay Index 

From the AIA observations at the time of the 
CME eruptions we deduce the initial directions of 
the CMEs in the low corona (directions shown in 
Fig[^. To understand why the initial directions of 
the two CMEs differ although they originate from 
the same AR, we calculate the decay index, i.e. 
the decay rate of the strength of the horizontal 
magnetic field with height for the entire computa¬ 
tional domain, defined as 


dlogBh 

dlogz 


(7) 


where = \JB"^ -I- B"^ and z is the radial dis¬ 
tance measured outwards from the solar surface in 
solar radii. The decay index is a measure of how 
the horizontal field strength of the magnetic field 
decays with increasing height. This quantity has 
been used to define a threshold for the triggering 
of ideal-MHD instabilities, i.e. the Kink (Tbrok 


& Kliem 2005) and Torus (Kliem & Torok 


2006) 


instabilities. For instance, for bipolar configura¬ 
tions and for n > 1.5, an MFR is Torus-unstable 
and thus can not be constrained by the overlying 
field and could erupt. In addition, the kink in¬ 
stability requires significant twist to be stored in 
the pre-eruptive magnetic configuration, ^ 1.6 
dTorbk fc Klie^[20^ . 


In Fig we show the distribution of the crit¬ 
ical decay index at different heights before and 
after the time of the CME eruptions. For the 
decay-index calculations we used the NLFFF ex¬ 
trapolations field discussed in the previous para¬ 
graphs. The locations at each height (shown with 
different color contours corresponding to different 
heights) form closed contours and the big con¬ 
centric “holes” denote the “weak spots” of the 
horizontal field, where the decay index is larger 
than the critical (or “super-critical”; here n > 1.5). 
This corresponds to areas with extremely weak 
horizontal component of the magnetic field (thus 
almost radial field lines), which is typically the 
case above the sunspots, e.g. at the leading 
sunspot in the west side of the AR. 


However, the magnetic conhguration may allow 
for similar ’holes’ away from the radial sunspot 
fields. For example, this is evident at the NE PIL 
(see family of black to blue contours), where a 
“cone” with decay index greater than the critical 
decay index is in close proximity to the MFR mag¬ 
netic structure along the PIL. Thus, a structure 
which is becoming Torus-unstable would preferen¬ 
tially erupt, and most probably “escape” through 
this super-critical decay-index “tunnel”. This can 
be seen in the three frames provided in Figure [TtI 
A movie covering the entire 2-day period of March 
6-7 is available in the online version of the paper 
(“movie5.mpg”). A comparison of the panels fo¬ 
cused on the NE PIL shows a migration and devel¬ 
opment of the supercritical decay-index “tunnel” 
stretched along the NE PIL (second panel). A 
comparison between the starting height of the “es¬ 
cape” tunnel (which before the eruption reaches all 
the way down to the surface), the cross-sectional 
radius of the NE MFR from NLFFF (^ 10 Mm) 
and the apparent initial height of the erupting 
structure from EUVI-B (Fig 10 ~ 30 Mm) yields 


consistent results supporting the necessary ingre¬ 
dients for an ideal instability to kick in. In ad¬ 
dition, the site of the initiation of CME2 corre¬ 
sponds also to low-lying supercritical decay-index 
regions. A comparison between Fig and the 
EUV observations shows good correspondence be¬ 
tween the observed direction for the CMEs, and 
the locations of both the supercritical decay in¬ 
dices and our presumed MFRs. The critical index 
analysis, therefore, justifies the discrepancy in the 
two CME directions and provides additional, indi¬ 
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rect support for the existence of MFRs since such 
structures are prone to plasma instabilities depen¬ 
dent on the critical index. 

5. Discussion and Conclusion 

We embarked on this study aiming to uncover 
the reasons behind this unusual twin CME occur¬ 
rence. The double CME events of March 7, 2012 
originated from a very dynamic and rapidly evolv¬ 
ing AR. Both CMEs were super-fast {vcme > 
2000 km/sec), launched within one hour of each 
other, and were associated with two major so¬ 
lar flares (X5.4 and XI.3). They, however, fol¬ 
lowed very different propagation paths. The 
Earth-directed CME, CME2, misidentified in pre¬ 
vious works, caused the second larger geomagnetic 
storm of Cycle 24. The source AR exhibited all 
the hallmarks of a “geomagnetically dangerous” 
AR: it had an anti-Hale d-spot magnetic configura¬ 
tion, multiple and compact PILs, strong shearing 
motions and non-potential shear. In other words, 
AR 11429 had plenty of magnetic energy and com¬ 
plexity to expel multiple, fast CMEs, as it clearly 
did during its disk passage. 

The continual flaring/CME activity of the 
source AR did not hinder our analysis, thanks 
to the detailed coronal coverage from SDO and 
STEREO. We were able to trace the origins of 
the double CMEs of March 7 to confined flar¬ 
ing events of various magnitudes on the previous 
day. For example, an M2.1 flare on March 6 at 
12:35 UT illuminated MFR-like structures to a 10 
MK temperature (in the AIA 131 A images). The 
position of the AR on the solar disk on March 6-7 
(E30°N18°) posed some challenge in our analy¬ 
sis. The AR structures were viewed from above. 
As discussed in the Introduction, the background 
emission (from low-lying structures) can impede 
the identihcation of MFR structures in such pro¬ 
jections. We focused on detailed NLFFF extrapo¬ 
lations using carefully treated HMI magnetograms 
over the full two-day period of interest. The ex¬ 
trapolations supported the presence of MFRs on 
March 6 as they showed MFR-like held lines— |J|- 
concentrations corresponding to oval-shaped held 
lines in vertical plane cuts (Figure [l^ — which 
are effectively weakly twisted structures with con¬ 
tinuous twist build up. 

The detailed NLFFF extrapolations allowed to 


estimate the poloidal hux of the MFRs be¬ 
fore the CME eruptions. For the NE MFR (as¬ 
sociated with CMEl) <l>p ^5.82x 10^° Mx and for 
the SW MFR (CME2) -3x lO^i Mx. The SW 
MFR, although shorter, showed higher spatial fre¬ 
quencies and thus more crossings in the calcula¬ 
tion plane per unit length. It was also rooted at 
stronger magnetic helds. The height for the cal¬ 
culation was determined to be up to 15 Mm since 
any further increase in height did not change the 
values considerably. This was consistent with our 
finding of low-lying MFR structures (Section |4.1[ ). 
Unfortunately, the MC httings for CME2 were 
highly uncertain so we were unable to derive an 
1 AU poloidal flux to compare with our origin es¬ 
timates. Comparing our fluxes with statistics of 
MCs observed at 1 AU we hnd good agree ment 
(e.g., ~ 5x 10^° - 2x 10^^ Mx, (Figure 8, Qiu 

et al.||2007[), $p ^ 8 x lO^o - 1.6x lO^i Mx (Table 


-LW i.UXN -LU J.VJ.W 

Mostl et al.|2009)). In particular, our value for 


1 , 

the poloidal flux seems to be at the center of the 
statistical distribution of poloidal fluxes from Fig. 
2 of Lynch et al. (2005). We surmise that a sig- 
nihcant amount of poloidal magnetic flux possibly 
existed in the MFR before its eruption. Overall, 
the combination of the solar observations and de¬ 
tailed NLFFF modeling enabled us to detect the 
MFR structures at least 12 hours before the erup¬ 
tions. We summarize our major hndings in the 
following list: 


1. The AR contains two major stressed mag¬ 
netic flux rope systems (MFRs), one in the 
NE and the other in the SW. We propose 
that the NE MER (giving rise to CMEl) 
formed during major confined flares, while 
the SW MER (giving rise to CME2) formed 
during confined sub-flaring episodes. 

2. We observe strong shearing motions. Oppo¬ 
sitely oriented magnetic elements move close 
to each other and give rise to sharp, com¬ 
pact polarity inversion lines (PILs); one PIL 
extends along the NE and another, shorter, 
PIL extends along the SW of the active re¬ 
gion (Fig|^& online movies “moviel.mpg”, 
“movie2.mpg”). These shearing motions last 
throughout our observing period and even¬ 
tually lead to a non-potential field topology. 

3. The strong shear leads to MFR formation 
via reconnection events during confined flar- 
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ing episodes, ranging from proper flares to 
sub-flares. The flare heating ’’illuminates” 
the newly created MFR field lines in hot AIA 
channels (e.g., standard-deviation maps of 

Fig[^. 

4. We find that a similar number of pixels in 
the vicinity of the PIL has both high twist 
(T„ > 0.5) and high shear {6 > 75°). This 
suggests weakly twisted non-potential struc¬ 
tures along the PIL. 


5. Further support that the NE MFR forms 
during the confined flare on March 6, is of¬ 
fered by the appearance of cool filamentary 
plasma along the presumed MFR axis dur¬ 
ing the cooling phase (~13:49:21). This in¬ 
dicates the existence of a coherent structure 
along the PIL, able to support cool plasma 
(in agreement with MFR observations of 


Patsourakos et al. (2013) that show cooling 
at the MFR center; see also “movieS.mpg”). 


6 . The pre-eruptive MFRs are low-lying and 
weakly-twisted. This is a challenge for their 
detection due to background emission in the 
lower corona. However we see them during 
flaring events when they do clearly stand out 
from the background. 


7. The low twist, seen in the EUV observations 
and produced in the NLFFF extrapolations, 
favors the ideal-MHD torus instability as the 
driver of the eruption. 


8 . The two CMEs originate from two PILs 
within the same active region. There is evi¬ 
dence for pre-existing MFRs over both PILs. 
The two PILs are two different flux systems 
and the formation of the two MFRs reflects 
that duality. The NE MFR is formed via a 
confined flare ~ 12 hours prior to its CME 
(e.g. similar structure in right column of 
Fig 12 hours apart; yellow dotted line). 
The SW MFR forms via small-scale flaring 
over many hours. The two CMEs continue 
to express their different nature even during 
the eruption where they follow paths deter¬ 
mined by the gradient of the magnetic field 
above each MFR. 


In conclusion, we were able to use detailed ob¬ 
servations from SDO and STEREO and sophisti¬ 


cated NLFFF extrapolations to decipher the ori¬ 
gins of the double CME event that occurred on 
March 7, 2012. They erupted from two differ¬ 
ent PIL systems, within the same active region. 
By tracing the coronal activity along those PILs 
backwards in time, we found strong evidence for 
the formation of magnetic flux ropes during con¬ 
fined flaring events of various magnitudes ranging 
from a proper flare (the M-class flare on March 6 
at around 12:35 UT) down to sub-flare events. 

Hot (~ 10 MK) elongated EUV structures ap¬ 
peared at those times. Their behavior was con¬ 
sistent with an MFR seen face-on considering 
past observations of MFR detections on the limb. 
NLFF extrapolations based on HMI observations 
over the full two-day period supported the EUV 
observations. They indicated the gradual forma¬ 
tion of low-lying, weakly twisted structures at the 
same locations where the hot EUV signatures were 
detected. Furthermore, the evolution of the large 
scale field structure around the AR exhibited two 
areas of sharp gradients in the magnetic field mag¬ 
nitude as a function of height. Such locations are 
considered the ’pathways’ through CMEs tend to 
erupt in the low corona since they provide the 
path of least resistance for the erupting system. 
True to this modeling result, both CMEs on March 
7, 2012 followed their corresponding ’pathway’, 
which explains their widely different propagation 
directions. In short, the two CMEs on March 
7, 2012 were the eruption of magnetic flux ropes 
formed 12 hours earlier during a confined flare 
episode (for CMEl) or over several hours via sub- 
flaring activity (for CME2). In other words, the 
two CMEs should be considered as separate events 
because they originate from different flux systems 
even though the two systems belong to a single 
AR. It raises the question on whether such com¬ 
plex (5-spot regions should be considered as single 
ARs or not. 

Finally, we note that the sheer complexity of 
the photospheric and coronal configuration of the 
studied AR presented a formidable challenge for 
the data analysis and NLFFF modeling. However, 
such complex, but rare, ARs, are the most dan¬ 
gerous in terms of space weather conditions and 
their study is therefore vitally important in under¬ 
standing the solar conditions that lead to extreme 
eruptive activity. 
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Fig. 1.— Full disk image from SDO/AIA at 131A identifying AR11429. The size of the box is 500"x400". 
The direction of the CMEs at the onset of the eruptions is shown with the green arrows. 
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Vector Magnetic Field of NOAA AR11429 



Fig. 2.— A sample CEA magnetic vector map from the HMI/SDO taken at the time indicated. The greyscale 
image is the normal component of Uphot saturated at ±2500 G. The horizontal photospheric magnetic field, 
B,* is shown with blue vectors. Note the alignment of the horizontal field along the PIT (observational 
manifestation of high PIL-shear). The two orange boxes shown enclose the NE and SW PIT respectively 
(boxes show FOVs used in Fig[T^. The green dotted window was used for the calculation of the magnetic 
flux (in Fig|^. 
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Magnetic Flux vs time for NOAA AR11429 
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Fig. 3.— Top: The evolution of the magnetic flux for AR 11429 between Mar 5 and 8, 2012 (Black diamonds: 
total unsigned, blue: positive, red: negative). The 1-day period oscillatory behavior is due to intrinsic 
instrumental problems of the HMI magnetograph becoming obvious in strong-field regions. Nevertheless, 
the overall trend is increasing in accordance with a continuously developing AR 11429. There are two 
distinctive flux emergence events showing their imprints in the flux profile - one on Mar 5, 2012 and a second 
one during/after the eruption of the twin CMEs (shown with arrows). Bottom panels: The time evolution 
of the photospheric vector field during March 5 and 6, 2012, i.e. before the eruptions of March 7, 2012. 
The blue vectors illustrate the horizontal component of the photospheric magnetic held. With red color we 
delineate the PIT (obtained by calculating the gradient of the component). Note that while the overall 
spatial distribution remains the same, significant changes occur along the PIT by means of shearing motions 
of individual flux elements (green dashed circles; mc^ notably in the NE side of the PIL - white “blobs”, 
following an emergence event on March 5, 2012 close to the NE negative polarity; location is shown with a 
red circle). Shearing occurs at the location of the SW negative sunspot with the elongated positive polarity 
on top of it. A movie (“moviel.mpg”) is available in the online version of the journal. 





















GOES X-Ray and AIA EUV Flux vs. Time 
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Fig. 4.— Soft X-ray flux versus time observed by the GOES satellite (black curve). The EUV normalized 
flux from AR 11429 for a few AIA channels is overplotted. The size of the box for the EUV flux calculation 
is 600" X 500" (corresponds to the box in Fig[^. Note the similarity in the response of the 131 A (cyan curve) 
with the 1-8 A from GOES which shows that AR 11429 is the primary contributor of X-ray flux during the 
observations. 
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Fig. 5.— Time evolution of 131A images from SDO/AIA. The thin super-imposed curve shows the location 
of the PIL at the photosphere (red shows location where the shear angle between the observed photospheric 
horizontal field and the potential field as a reference, is greater than 60°). The PIL was calculated at the 
CEA Bz frames of fig. [^and it has been transformed into helioprojective to be used with the 131 A images. 
The majority of the transient activity resides at the PIL, and due to the position angle (~ E30°N18°), the 
transient bright structures are close to the photospheric surface. The MFRs are delineated with a dotted 
curve; yellow for the NE MFR and pink for the SW. The NE and SW MFR are offset to the PIL, which 
signifies potentially higher lying structures. 
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Fig. 6.— The evolution of the AR in AIA/131A with HMI/BLOS composite (gold/positive, navy 
blue/negative). The locations of the NE MFR and SW MFR are marked by arrows in the middle top 
and bottom panels. During the cooling phase of the failed eruption 13:48 UT), an absorption feature 
appears along the presumed NE MFR axis. The absorption feature is offset from the PIT, which may suggest 
that it is a high lying formation (due to the geometrical projection from our viewing angle). This abso rption 
feature is also seen in 195 A from STEREO B’s viewpoint and it is indeed high-lying (see Fig 10). The 
associated movie to this figure (“movie2.mpg”) can be found in the online version of the paper. 
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Fig. 7.— 211A base difference images around the time of the M2.1 flare of March 6, 2012. The absence of 
spreading dimmings areas in the vicinity of the AR suggest that this was an confined flare. 
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Fig. 8.— Top six panels show the MFR candidate structure in the six coronal passbands (131 A, 94 A, 
335 A, 211 A, 193 A and 171 A) observed one hour after the confined M2.1 flare of March 6, 2012. In the 
bottom six panels (colored in green) we present DEM images in the displayed temperature ranges. The time 
of each temperature frame corresponds to the mean time of the coronal passbands (maximum time difference 
between passband images is 10 sec). The MFR candidate is seen above 2MK (dotted line). 
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Fig. 9.— Left: Slit time plot for EUVI B 195 A (top), AIA 193 A (middle) and 131 A (bottom). The EUVI 
B slit plot’s time cadence is 300 sec and the AIA’s is 12 sec. The vertical axis is in arcseconds along the 
slit. Right: Selected frames showing the location of the slits for EUVI B and AIA images. The slit was 
selected to intercept the image along the direction of fastest expanding motions. The frames in the last 
column are laplacian-filtered to enhance the details. A movie is available in the online version of the paper 
(“movie3.mpg”). 
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Fig. 10.— The final stages of the failed eruption of March 6, 2012 12:30 UT as seen by STEREO B/EUVI 
at 195 A. Note the absorption feature along the presumed axis of the NE MER shown with red arrows. This 
is evidence for a magnetic structure able to hold cool material as expected for the axis of the MER. The 
height of the absorption feature is ^30 Mm. The formation of the absorption feature observed in 131A can 
be seen in the online movie (“movie3.mpg”). 
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Fig. 11.— NLFFF extrapolation of March 6, 2012 23:48 UT. The greyscale map is the component of 
the photospheric boundary at that time. The colored tubes represent the extrapolated magnetic field in the 
domain. The teal-colored tubes correspond to the field lines along the PIL. The magnetic field lines are 
highly sheared which indicates the existence of strong electric currents in the vicinity of the PIL. These field 
lines are rooted in randomly sampled points within areas of |a| = 5 x at the surface. A movie 

is available in the online version of the paper (“movie4.avi”). Note the existence of two main “chains” of 
non-potential field lines - a short one, above the negative sunspot in the South and a longer one following 
the NE PIL. These correspond to the locations of the brightenings seen in 131A and also to the locations of 
initiation of the two CMEs. The green tubes represent the overlying (nearly potential) field lines sampled 
randomly in the FOV. 
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Fig. 12.— Top row: 131A images showing the evolution in coronal plasma emission of AR11429 at 10 MK. 
Middle row: Volumetric visualization of the magnitude of strong (normalized to peak current density) 
electrical currents calculated via Ampere’s law in the NLFFF volumes. The FOV roughly corresponds to 
that of the top row; the volume structure is seen from above. The volumetric currents are represented 
in an “optically-thin” fashion and a teal color table has been used for the visualization. These strong 
currents reside close to the surface as also seen in the vertical plane cuts (see bottom row). Bottom row: 
Vertical cuts of the 3-D volume |J| at the position indicated in the middle row (black line). Color contours 
show strong |J| with red and weak with blue. The magnetic field cut by the vertical plane is overplotted 
(note the strong curl in regions of high |J|). The greyscale maps correspond to the photospheric boundary. 
Transient brightenings in the 10MK (131 A) channel seem to be in good agreement with |J|-volume channels 
inferred from the NLFFF. The NE PIL appears fragmented at earlier times (white arrows) and progressively 
straightens, lengthens and becomes an integral structure as signified by the transient brightenings. Several 
hours prior to eruption, two main hot |J|-channels are visible, in agreement with 131A . 
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Fig. 13.— Helicity-related parameters from the extrapolation of March 6, 2012 23:46 UT (i.e., the same 
shown in Fig [TI|). Top four panels: NE of fig[^: (a) Bz map, (b) force-free a map, (c) map of field line 
lengths, L, and, (d) the unsigned twist number |T„|. PIL is shown with a white curve. Bottom four panels: 
Same parameters but for the SW box of Fig Only a small fraction of pixels in the vicinity of the PIL 
exhibits significant twist (|Tji| > 0.5). 
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Fig. 14.— The 2-day evolution of the average force-free a (top), average twist number (middle) and the 
fraction of flux with twist higher than a twist threshold (bottom). Note that the y-axes for top and middle 
panels are reversed. Only a fraction of the vicinity of the PIL exhibits significant twist. The twisted flux 
increases after the two CMEs (vertical line), which suggests continuing flux emergence past the time of the 
CME eruptions. The unsigned magnetic flux is shown in the middle panel (red). Note the fast increase of 
magnetic flux right after the eruptions. Recall that our results are at 1-hour cadence, which compares to the 
temporal lag between the onsets of the two CMEs. 
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Fig. 15. — Composites of 13lA observations and force-free a-maps calculated at the photosphere. Left panel 
shows the map at March 6, 2012, 10:59:56 UT and on the right at March 6, 2012, 23:47:56 UT. Blue/red color 
corresponds to positive/negative a. The teal-colored overlaid map corresponds to the ^2-hour accumulated 
brightenings (up to 2-hours before the aforementioned times) in the 131 A AIA passband by calculating the 
standard deviation of each pixel in the 131 A image time series. Note that the a-maps are converted back 
to helioprojective coordinates to match the AIA standard deviation maps. 


30 








00:00 08:00 16:00 00:00 08:00 16:00 
Start Time (05-Mar-12 23:29:56) 



0 66 131 00:00 08:00 16:00 00:00 08:00 16:00 

Start Time (05-Mar-12 23:59:56) 


Fig. 16.— Left: The distribution of the non-potential shear angle 9, i.e. the angle between the horizontal 
component of the observed photospheric field, and the held from the potential extrapolation, Bpot, 

at an hourly step, for 2 days (each column of the left panel of this figure corresponds to the theta angle 
histogram at the corresponding time). The mean shear builds up ~ 4 hours prior to the X5.4 flare of 00:02 
UT March 7, 2012 and the 2cr-range of shear shifts with the mean shear. The fact that the shear is relatively 
strong before and after the CMEs, may suggest that part of it is involved in the MFR formation. Right: 
Total number of pixels in the PIL band above a 0-threshold. For 9 >75° the increase of the number of highly 
sheared PIL pixels seems to occur ~ 8 hours prior to the eruption. 
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Fig. 17.— Decay index with height contour plot overlaid on photospheric CEA Bz maps (color contours 
correspond to decay index n=1.5 at each height). The decay index profile is calculated for each z-column 
of the horizontal magnetic field from the NLFFF cubes. Top panel is at the beginning of the observations, 
middle is a few minutes prior to CMEl and bottom right before CME2. Areas within color-contours contain 
“super-critical” decay indices at each height, n > 1.5. Note the activity in the low-lying (darker colors) 
contours of the NE PIL. In just one day (top frame to middle) an elongated and continuous stacking of 
“super-critical” decay-index areas is being formed, giving birth to a “super-critical” decay-index “tunnel” 
(see succession of closed contours at the location pointed by the red arrows). The succession of these super¬ 
critical decay-index “holes” signifies a least-resistance path, or the most probable path for the CMEl initial 
trajectory (i.e. towards the East). A movie is available in the online version of the paper (“movieS.mpg”). 
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Table 1 

MFR Length, Cross-section, Core Height and mean |B| before the CME eruptions* 


MFR 

Length (Mm) 

Cross-section (Mm)l 

Core height (Mm)^ Mean |B| (G) 

NE MFR 

80 

10 

3 

650 

SW MFR 

40 

6 

3 

500 


* Determined from the NLFFF extrapolation of March 6, 2012, 23:48 UT 

Cross-section inferred from the 20% of peak |J| at the slice cut along the longest 
direction. 

tHeight inferred from peak value of |J|. 
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A. NLFFF Extrapolation 


The boundaries were binned twofold (2x2 binning, i.e. down-sampled to 282x 195 pixels) to accommodate 
our computer resources. A preprocessing of these rebinned photospheric boundaries followed, to make them 
compatible with the force-free condition (as in [Wiegelmann et ah 2006). The size of the computational 
domain is 282x195x118 pixels and a potential magnetic field extrapolation was carried for each boundary 
of our time-series (using the component only) to serve as the initial field to be optimized into a 

NLFFF field. The potential extrapolations were calculated via a standard Green function integrator for the 
magnetic scalar potential, Then the potential datacubes, derived as 'Rpotn = are fed into the 

NLFFF optimization routine along with the (preprocessed) photospheric boundaries 
The boundary layer at each of these faces was chosen to be 16 pixels thick. This resulted in a “physical” (i.e. 
numerical boundary-effect-free) sub-volume with a size of 250x163x106 pixels (or 180xll7x76Mm^) for 
each component of the NLFF magnetic field. This is the size of the final products of the NLFFF extrapolation 
used in our analysis. 


B. Twist Number Calculation 

Since points-of-avoidance or holes on the length-map signal overlapping field-line endpoints in neighboring 
pixels, we modified the generation of length-maps in the following way. We increased the sampling step for 
the seedpoints ten-fold using bi-linear interpolation and stored the locations of footpoints with the value for 
the length in a 10x larger length-map. This allowed for a 10-fold increase of the “resolving power” of the 
image grid for discriminating the locations of the end-points of field lines. This large length-map is therefore 
sparse and in order to fill the gaps we performed a forward integration from the locations of the holes at the 
positive polarity (end-points from the initial integration). Then the (now negative) end-points are stored 
in the same map and a Delaunay triangulation is performed. The latter is necessary in order to fill the 
remaining gaps within the negative and positive polarities via bi-linear interpolation. Then, multiplying 
this resulting length-map with its associated (i.e. 10 x larger via bi-linear interpolation) d-map, over 47r, 
yields the twist number, Tk, for the NLFFF extrapolated field lines in the physical domain. In this way we 
can derive smooth and continuous r„ maps for each of the NLFFF datacubes and now can study the time 
evolution of field-line twist as a way to monitor topological changes in the NLFFF. 
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